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Abstract 



. The asymptotic analysis of a linear high-field Wigner-BGK equation 

is developped by a modified Chapman-Enskog procedure. By an expan- 
sion of the unknown Wigner function in powers of the Knudsen number 
. e, evolution equations are derived for the terms of zeroth and first order 

in e. In particular, it is obtained a quantum drift-diffusion equation for 
the position density, which is corrected by field-dependent terms of order 
e. Well-posedness and regularity of the approximate problems are estab- 
lished, and it is proved that the difference between exact and asymptotic 
solutions is of order e 2 , uniformly in time and for arbitrary initial data. 



Key words: Asymptotic analysis, quantum drift-diffusion model, Wigner 
equation, open quantum systems, singularly perturbed parabolic equations. 



"*-! 1 Introduction 



Quantum mechanics has recently proved an essential tool for modeling the new 
generation of nanodevices [30]. However, the adoption of quantum models re- 
quires a delicate compromise with quantum statistics principles. Hamiltonian 
dynamics is described at the quantum level, either in terms of wave-functions 
(via Schrodinger-Poisson-systems), or of density- matrix operators (via von Neu- 
mann equation). For different reasons, both formulations are not suitable for 
simulations: precisely, the wave-function approach can not be extended to pic- 
ture dissipative dynamics of open quantum systems, while the density matrix 
approach is not appropriate to describe finite position domains, due to its non- 
local character. For the same reasons, it is instead convenient to employ (Wigner) 
quasi-distribution functions [6l 133]. Nevertheless, a phase-space description of a 
multi-dimensional dynamics presents well-known computational drawbacks. On 
the other hand, quantum hydrodynamic models seem to be a promising tool 
both from the numerical and the analytical point of view [22l [23] . Similarly, in 
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semi-classical semiconductor theory, the interest of modelists has shifted from 
Boltzmann equation to hydrodynamic systems, and they have been widely stud- 
ied both for a physical validation and from an analytical and numerical point of 
view (cf. [1] and the references therein). A rigorous derivation of quantum hydro- 
dynamic models from more fundamental ones, in either Schrodinger or Wigner 
formulation, is an open and analytically demanding problem [HI |2"TI |2"3"] . This is 
the motivation of the present paper. 

The preliminary step for passing from the kinetic picture to a macroscopic one 
consists in including dissipative mechanisms in the evolution model, for exam- 
ple, the interaction of the quantum system with the environment. In the weak 
coupling limit, a Markovian dynamics can still be adopted, and the description 
via an (operatorial) evolution equation in Lindblad form is considered quantum- 
physically correct [26]. From this class of evolution equations, kinetic models of 
open systems can be derived via Wigner transform. In Section 2 we shall briefly 
review the most popular Wigner models of irreversible dynamics. 
In this paper, we consider the case of an open quantum system in a high-field 
regime, more precisely, of an electron ensemble subject to an external poten- 
tial, whose effect is comparable with the interaction with the ion crystal. In- 
cluding high-field effects has great relevance in semiconductor simulation. A 
macroscopic model of this evolution is expected to contain field-dependent trans- 
port parameters, that are tipically deduced via fitting procedures. We refer the 
reader to [13] for an updated review of derivations of semi-classical high-field 
drift-diffusion models by diverse limit procedures: in particular, in [13], from an 
energy-transport model, are obtained explicit field-dependent mobilities. On the 
contrary, in [9] , a high-field drift-diffusion model with non-explicit field-dependent 
coefficients is derived, as the limit of a Spherical Harmonics Expansion of semi- 
classical Boltzmann equation. 

We present a rigorous derivation of a Quantum Drift-Diffusion (QDD) equation 
with explicit field-dependent mobility and diffusion coefficient. We shall start 
from the Wigner equation with an additional linear BGK term, modeling the 
interaction with the environment, and then adapt the equation to the high-field 
case, by rescaling it in terms of the Knudsen number e. Thus, our contribu- 
tion is the quantum counterpart of [52] • We recall that, in [20], the starting 
point is the Wigner-BGK equation as well, but collisions are considered to be 
the strongest mechanism during the evolution (moderately high-field regime), 
and the relaxation term is derived via a Chapman- Enskog procedure. In our 
case, the additional relaxation term is instead an 0(H 2 )— approximation of the 
Wigner-transformed relaxation term in operatorial formulation (cf. Section 2). 
Moreover, we perform an asymptotic expansion of the unknown Wigner function 
in terms of e, according to a modified Chapman-Enskog procedure introduced in 
[3 lj . This method has been applied to many kinetic models and constitutes a 
valuable tool for a rigorous asymptotic derivation of macroscopic models (cf. Sec- 
tion 5). We substitute the Wigner unknown in the originary evolution problem, 
with the expansion of order e 2 , and we get an approximated problem: in par- 
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ticular, an equation with unknown the electron position-density. This equation 
is precisely the QDD equation corrected by the (9(ft 2 )-Bohmian term of order e, 
and by field-dependent terms, of order e as well. These terms contain the same 
field- dependent coefficients obtained in the semi-classical case [T3| [32]. 
The well-posedness of the 0(e 2 )-approximated problem is discussed in Sections 
7 and 8, and finally, in Section 9, we prove that the difference between the so- 
lutions of the originary and of the approximated evolution problems is also of 
order e 2 . In conclusion, with the present analysis we obtain a QDD equation 
with field-dependent mobility and diffusion coefficients and we prove rigorously 
that, up to a certain degree of accuracy, it constitutes a model of quantum trans- 
port in the high-field case. From the analytical point of view, this equation is a 
second-order parabolic PDE with non-homogeneous coefficients. In particular, 
it belongs to the class of singularly perturbed equations; accordingly, the well- 
posedness result, together with the regularity estimates derived in Section 8, are 
complementary to the discussion in [8] about the same class of equations with 
constant coefficients. A counterpart of our analysis is the well-posedness study 
of the quantum drift-diffusion equation, in the fourth-order formulation obtained 
via a "classical-equilibrium" approximation [23] . We remark that the asymptotic 
procedure used here presents analogies with the Chapman- Enskog one in kinetic 
theory; nevertheless, it is well-known that the latter does not deal with the "initial 
layer" problem, namely, the instants close to the initial one are excluded from the 
analysis, due to the rapid changes of the solution [32] ■ In the present approach 
instead, the initial layer problem is solved at once. 



2 Wigner-BGK equations 



Let us consider a quantum system with d degrees of freedom, evolving under 
the effect of an external potential V = V(x), x G IR d . The Wigner equation 
with unknown the quasi-distribution function w = w(x,v,t), (x,v) G IR M , t > 0, 
provides a kinetic description of the evolution of the system. It reads 



dw 
~dt 



+ v-V x w-Q[V]w = 0, (x,v)gIR , t>0 



(2.1; 



with the pseudo-differential operator Q[V] defined by 



(Q[V]w)(x,v,t) 



(2yr) d 
% 



EVE 



E d 



6V(x, rj)w(x, v', t)e l{v - v ' yv dv' dr) 
5V(x,r])Fw(x,ri,t)e lv - r >dri, (2.2) 



where 



5V(x,rj) := - 



V [ x- — 

2m 



and Tf{rf) = [T v ^ v f\ (rj) denotes the Fourier transform of w from v to r]. In 
the Fourier-transformed space IR^ x IR^ the operator 0[V] is the multiplication 
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operator by the function i 5V; in symbols, 

T (Q[V}w) (x, ri) = i SV{x, T))Fw(x, rj) . (2.3) 

Eq. (12. ip corresponds via Wigner-transform to the von Neumann equation de- 
scribing the conservative dynamics of an isolated quantum system [33]. Suc- 
cessive modifications of the Wigner model have been proposed to picture an 
irreversible interaction of the system with the environment. In [T7] a scattering 
term is derived by a weak-coupling limit; however, due to its non- locality, it is 
not suitable for simulations and for mathematical analysis. A second possibility 
is an additional diffusive term, as in the quantum counterpart of Fokker-Planck 
(FP) equation of classical kinetic theory [IT] (cf. [12] for the latest derivation and 
[HE] for the latest well-posedness results). Unlike the Wigner equation with the 
scattering term, the quantum FP equation is the Wigner-transformed version of a 
Markovian master equation in Lindblad form, namely, it is the kinetic version of a 
quantum-physically correct model [7] . The shape of the drift-diffusion equations 
corresponding to the low-field, respectively high-field, scaling of the classical, re- 
spectively quantum, FP equations are presented in [3J. Another possibility is to 
insert a BGK operator, either linear or non-linear, like in [TD] , meaning that after 
a time \ jv the system relaxes to a prescribed state w cq , namely, 

+ ^xW - e[V)w = -v(w - w cq ) , (x, v) e JR 2d , t > 0. (2.4) 

In the recent literature [TH [TU [20l [23], diverse relaxation-time states w eq have 
been proposed. 

The standard picture is that the system converges to a state of thermodynamical 
equilibrium with the surrounding environment at temperature T. The opera- 
tor that individuates the statistical equilibrium state at (constant) temperature 
T = l/k/3 (k is the Boltzmann constant) is e -1311 , H being the energy operator 
associated to the system. The von Neumann equation modified by a relaxation- 
time term containing e _/3H is in Lindblad form [2]. Accordingly, a Wigner-BGK 
model, being the Wigner-transformed version of that equation (i.e., containing 
the Wigner-transformed of e _/3H as relaxation-time state), formally belongs to 
the class of quantum-physically correct kinetic models. In his pioneer article 
[33] , E. Wigner applies an expansion in terms of h to the Wigner function corre- 
sponding to the operator e _/3H , and obtains the classical equilibrium distribution 
function on the phase space with correction of non-odd order in h: 
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— AV + — |W| 2 + /3 > v r v s 



m m ' dx r dx s 

r,s=l 



+ 0(h 4 ) , (2.5) 



where S(x, v) := mv 2 /2 + V(x) is the total energy of the system. Let us call w eq 
its local (in time and space) version, defined by 

w cq (x,v,t) := C(x,t) w w (x, v) , 
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with C to be chosen. By assuming 



w eq (x,v,t)dv = Jw(x,v,t)dv =: n[w](x,t) = n(x,t) 
and since, by the direct computation, 



(2.6) 



^w(x,,)^=(^) d/ >ll + n 2 ^ 



12m 



-AV + ^\VV\ 2 



the local Wigner thermal equilibrium function w eq equals 



Weq (x,v,t)=n(x,t) (^y 2 e-" mv2/2 



i + h 2 ^- 

24 



1 d 

—AV + py « r i 



r,s=l 



«9 2 V 

decree q 



+ 0(h 4 



(2.7) 



In ( 12 .7j) can be recognized the classical (normalized) Maxwellian 

:= 



f3 m \ d/2 e -(3mv 2 /2 



2tt / 



parametrized by the density n and the constant temperature l/k/3, with an addi- 
tional correction term of order h 2 . We shall consider the expression (12.71) as the 
(9( ^-approximation of the Wigner function associated to the state to which the 
quantum system shall approach. 



In [15] is presented an alternative strategy to individuate the relaxation-time 
state, that is the extension to the quantum case of Levermore's one for classical 
kinetic equations ([23], cf. [T] for semi-classical equations). It consists in tack- 
ling a constrained minimization problem for the relative entropy of the quantum 
system under consideration, with respect to the environment. In the quantum 
case the procedure is performed at the operatorial level, due to the non-local 
definition of the entropy, in terms of the operators describing the states of the 
quantum system. However, the constraints for the minimization procedure are 
considered at the kinetic level. Thus, the Wigner transform W is used inten- 
sively to pass from the operatorial formulation to the kinetic one, any time it is 
required by the procedure. Due to that, the expression of the minimizer of the 
entropy formally derived in [H] is non-explicit. Nevertheless, in [T3], is formally 
proved that W{exp W" 1 /} = exp/ + 0(h 2 ) with / defined on the phase-space. 
Accordingly, the (formal) minimizer reads 



w A (x,v,t) 



^A-/3mv 2 /2) 



1 + 



1 . . p 1V7 , l2 p A 
— AA + — \VA 2 + - > 
m 3m 3 ^— ' 



v r v< 



r,s=l 



d 2 A 



+ 0(h 4 )} (2.9) 
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with A = A(x, t) Lagrange multiplier used for the constrained minimization pro- 
cedure, i.e. 

J ' WA(x,v,t) dv = n(x,t) . 

By comparison of the expression (12.91) with (12.71) . it can be easily seen that they 
coincide if one identifies the Lagrange multiplier A with —j3V. In [23] it is indeed 
proved that A = -f3V + 0{h 2 ) holds. 

Remark 2.1 It is crucial to recall that the correspondence via Wigner-transform 
of the operatorial and the kinetic formulations is merely formal, unless certain 
assumptions are posed both on the Wigner functions and on the operators [27] . 
On this point depends the analytical difficulty in stating rigorously the well- 
posedness of the strategy of derivation in |14j . For the same reason, the analysis 
of Wigner equations is set in the Hilbert space L 2 , since the necessary condition 
for the rigorous correspondence is satisfied [27] (cf. [5], e.g.). □ 

As a consequence of the previous discussion, in the present article we shall adopt 
the Wigner-BGK equation f!2.4[) containing (12. 7p on the right-hand side as the 
model of the open quantum system evolution. In particular, we remark that we 
shall consider the operator on the right-hand side as an (9(/i 2 )-approximation, 
in the kinetic framework, of the dissipative dynamics induced by the interaction 
with the environment. 



3 The high-field Wigner-BGK equation 

Our aim is to describe an open quantum system subject to a strong external 
potential; in particular, the action of the potential is to be considered comparable 
with the interaction with the environment. In order to adapt the Wigner-BGK 
equation (12. 4p to this specific case, we rewrite it by using dimensionless variables 
and, for this purpose, we introduce the time-scales of the action of the external 
potential and of the interaction with the environment. Let us call ty the potential 
characteristic time and tc the mean free time between interactions of the system 
with the background. Then we introduce x' = x/xo, v' = v/vo, t' = t/t , with 
xo, vo, t characteristic quantities, and we call w' = w(x', v', t') the rescaled Wigner 
function (observe that we can indeed neglect to rescale the Wigner function). 
Thus, we obtain 

— ^■w + vV x w-—e[V]w = -—u(w-wJ), t>0, (x,v)eH 2d , 
v t at v t v v t c 

where we have omitted the prime everywhere. If we introduce the relation xq = 
voto, we obtain 

d t t 

— w + v ■ V x w- — e[V]w = — -u(w-w cq ), t>0, (x,v)EJR 2d . 

Ot ty tc 



6 



In the following we assume that the times ty and tc are comparable, in the sense 

ty tc 



e, (3.1) 
to to 

where e := 1/xq is the Knudsen number, since I := Vote is the characteristic 
length corresponding to the classical mean free path. This corresponds to say 
that the external potential and the interactions coexist during the evolution. In 
particular, e ps corresponds to an evolution in which the effect of the interactions 
is dominant on the transport (tc « to or equivalently / <K Xo). However, at 
this time the action of the external potential has the same strength, due to the 
assumption (13. 1 p [ty « t ). In fact, the resulting equation is 

e— + ev • V x w - e[V] w = -u(w-w eq ), t>Q, {x,v)eJR 2d . (3.2) 

We recall that it is the quantum counterpart of the one studied by F. Poupaud 
in [32]. 

Now we put (13.21) in abstract form. As motivated in Remark l2.ll a suitable setting 
for problems in Wigner formulation is the Hilbert space L 2 (IR M ). However, in 
order to give a rigorous sense to the expression 

n(x) := w(x,v)dv, Vx G IR d , (3.3) 



which enters the equation via the definition (12.71) of w cq , we introduce the subspace 
X k := L 2 (1R 2d , (1 + \v\ 2k ) dxdv; IR), with k G IN, endowed with the norm 

II^Hl. = / \ w (x,v)\ 2 (l + \v\ 2k )dxdv. 

Let us call X v k the Hilbert space L 2 (H d , (1 + \v\ 2k ) dv- IR) and if™ the Sobolev 
space H™ ®X^. The weight k has to be chosen according to the space dimension 
d: we call <i-admissible 

k G IN such that 2k > d . 
The definition (13. 3p is well-posed for all w G Xk with <i-admissible k, since 

/ w(x,v)dv <C(d,k)( f \w(x,v)\ 2 (l + \v\ h ) 2 dv) , Vxe IR d 

by Holder inequality (cf. [29]). We define the streaming operator S by 

Sw = -v V x w , D(S) = {w G X k \Sw G X k } , 
and the operators 

Aw:=Q[V]w, Cw:= -{vw -flio) , Vw G Xk (3.4) 
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with the operator Q defined by 



Qw(x,v) := uF{v) { 1 + h 2 ^— 



24 



1 d 
-AV + /3 V 
m z — ' 



v r v< 



r,s=l 



d 2 V 

doc f oc a 



w(x, v') dv' . 



The function F(v) is the normalized Maxwellian, given by (I2.8p . Observe that 
we substitute the function w eq defined in (I2.7P with the operator fho, that differs 
from w cq by terms of order h 4 . Let us call F^ the 0(h 2 ) — coefficient in the above 
definition of O 



d 2 V 

QXr> OC c 



F(v) 



F^{x,v)=F^[V]{x,v) = --AV + P V v r i 

24 m 

r,s=l 

such that 

flw{x,v) =vn[w](x) [F(v) + h 2 F {2 \x,v)] . 

Observe that such expression for Qw can be seen as an 0(h 2 )— correction to the 
classical product n(x)F(v). 

In conclusion, we write Eq. (I3.2I) in the abstract form 

dw 



dt 



e Sw + Aw + Cw, 



(3.5) 



lim f ^ + \\w(t) - w \\ Xk = 



where Wq is the initial condition. 

In next Lemma we specify under which assumptions the abstract definition 
(I3.4p of the operator A + C is well-posed. 

Lemma 3.1 If V £ H% with d- admissible k and AV G , then the operator 
A + C is well-defined from Xk into itself, and is bounded by 

\\A + C\\ B(Xk) < C(d, k) [\\V\\ H , + v\\AV\\ L ~ \\F\\ Xh2 + u\\F\\ Xjl + v 

Moreover, A + C is well-defined from XJ! into itself, and is bounded by 

\\A + C\\ B{xl) < C(d, k) (\\V\\ m + u\AV{x)\\\F\\ xy _ + v\\F\\ x , + v 



fc + 2 



(3.6) 



Proof. Here and in the following we indicate with C non necessarily equal 
constants. 

The arguments are similar to those in |28j . so we just give a sketch of the proof. 
First of all, by the product shape of the pseudo-differential operator in Fourier 
variables (cf. (12.31) ). for all w G X^, it holds 

2 



|6[V>| 



C\\5VJ 7 w\\j 2 +C 



T 2 

^x ,77 



< 2C\\V 



\w\ 



+ C 



d ak 



T 2 



since 2k > d guarantees that H% ^ L£°. Here, the constant C is due to the 
Fourier transform. Then, by applying the product-formula rule and using Sobolev 
embeddings for the functions V and Tw G L 2 X ®H^, it follows that ||Q[V]||b(Xj,) < 
C|| V||//it. Moreover, for all w G Xk with 2k > d, 



\n w \\ 2 Xk <v 



\2k\ 



M- 



p 4 h 4 

24 2 m 2 



\AV\ 2 (x) \F(v)\ 



1 2fc> 



2A 2 m 4 



<vC{l + \\AV\\ 2 )\\F 



r,s=l 

\xv\ 



d 2 V(x) „. 
v r v, q — — - — F{v 



\w\ 



x k +vC\\AV\\ 



2 



w(x, v') dv' 



dx dv 



w(x, v') dv' 



SR" 



fe+2 1 



\w\ 



dx dv 



(3.7) 



since F G Xk , V/c . Then the estimate of is straightforward. The esti- 

mate (I3.6p in XI can be proved analogously. ■ 

We remark that the existence and uniqueness of a solution in X^ of the initial 
value system (I3.5P for any e > can be stated under the assumptions of Lemma 
13.11 by using arguments of semigroup theory, analogously to [29J . 



4 Well-posedness of the problem with e = 

The aim of this paper is to perform an asymptotic analysis of the system f l3.5|) . 
by using a Chapman-Enskog type procedure. The first step of the analysis is to 
solve Eq. (13. 5p with e = 0. This corresponds to individuate the Wigner function 
describing the state of the system in case the interaction of the environment and 
the action of the potential are dominant with respect to the transport. We remark 
that the function w eq defined by (12.71) describes the state to which the system 
relaxes under the sole interaction with the environment. 

We consider the equation (A + C)w = in the space X^. the variable x can be 
considered as a parameter in the analysis, thus we shall study (A + C)w = in 
the space XI for any fixed x G !R d . However, with an abuse of language, we shall 
indicate the operators with the same letters also when x is fixed. We can state 
the following proposition. 

Proposition 4.1 IfVE H% with k = max{2, k}, and d-admissible k, then for a 
fixed x G H d 

ker(^ + C) = {cM(v),ce H} C X v k , (4.1) 



with 



M(x,v) 



v — i5V(x, 7]) 



(3h 2 

24m 2 



r,s=l 
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for any fixed x. Moreover, for all h £ (A + C)w = h has a solution if and 
only if 

/ h(v) dv = . (4.3) 

Remark 4.1 It can be immediately deduced by the characterization (14. ip that 
the solution of the equation (A + C)w = in X k is unique, except for a factor of 
the sole x. □ 

Proof. By definition, 

ker(^ + C) := {w £ X v k \ [y - B[V])w = Q w}. 

For all h £ X k , the Fourier-transformed version of [y — 0[V])u> = h reads [y — 
iSV)J r w = Th. Thus, 

w(v) = {v- Q[V}Y l h{v) := T' 1 ( u ™gy {v) ) ( v ) ( 4 - 4 ) 

is the unique solution; equivalently, the operator (y — 0[V]) is invertible in XI 
with bounded inverse, defined by (I4.4p . Precisely, 

IMIi r IK- - e[v]rhf x , = c j-^0L_ ir, + c±J 

s~i dp 
r=l J 

then, by applying product formula, it can be checked that, if 2k > d, 

\\(»-®[V})-% ( x V <C(l + \\V\\ Hk J. (4.5) 

Then, 

w £ ker(^l + C) w = (u - QiV^ttw 

w = un[w}(u - efV])- 1 ^ + h 2 F {2) [V}) w = n[w]M(v) , (4.6) 

by definition of the operators iy — OfV"])" 1 and Q . From this follows the charac- 
terization (14.11) of ker(^4 + C), with the function M defined by 

M(x,v) := v{y - BSy\Y x {F{y) + h 2 F {2) [V}(x, v)) \/v £ IR d , 

with the fixed x £ IRA Since F+h 2 F^[V] £ X k for all k £ IN, provided AV £ L 2 X ; 
then, due to the assumption on V and to (14.51) . M £ X k if 2k > d. For all h £ X k , 
solving („4 + C)w = h is equivalent to (J — (u - e[V})~ 1 n)w = -(u -Q\V\y l h. 
Moreover, by the equivalence (14.61) . ker(^4 + C) — ker(X — [y — Q\y\)~ x VL ). Since 



d k Fh{r)) 
drjr v — i5V{rj) 



drj 



d k Th{y){y + iSV{r])) 

CLl I 



drfc v 2 + (8V) 2 (r]) 
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ker(.4+C) ^ {0} , the operator I—(u — Q[V]) 1 fl is not injective. If the operator 
[y — 6[V]) _1 f2 is compact, by the Fredholm alternative, this is equivalent to 
R((u - e jv])" 1 !] ) ^ XI The equation (I - (u - OfV])" 1 ^ )w = M has indeed 
no solution, since 



(1 - (u -Q[V]y 1 Q)w( y v)dv = 0, WweXl 



k i 



JK" 



(by the definition of the operator [y — 0[V]) while, instead, J M(v)dv = 

J F(v)dv = 1. Analogously, for all it e ker(.4 + C), (I — (v — Q[V]y 1 n)w = u 
has no solution. 

In conclusion, if we show that [y — OfV])" 1 ^ is a compact operator, then we 
can conclude by the Fredholm alternative that (X — (y — 0[V]) -1 f2 )w — h has a 
solution iff J h(v) dv = 0. Analogously to Lemma 1 in [32], it can be constructed 
by Rellich-Kondrachov theorem, a sequence of bounded finite rank operators 
converging to [y — G[V])~ 1 fi . Thus the thesis follows. ■ 
Finally, let us compute the first and the second moments of the function M: 

Lemma 4.1 Let V G H^ +2 with d- admissible k. Then, the function M defined 
by (14. 2 p satisfies 

(A + C)w = & w = n[w}M with J M{x,v)dv = l, (4.7) 



vM(x, v) dv 



vm 



■W(i) 



(4- 



r 12 Bh 2 

/ v ® vM(x, v) dv = — + ^— r W ® W + i^—r V <g> W . (4.9) 

J pm v^m 1 12m z 

Proof. ( 14. 7D follows by (14.6p . Moreover, since V e #£ +2 with 2fc > d, then the 
function M 



M(x, v) = v. 



T _ x ( r(F + h*FV) \ 
\ v — i5V(x, rf) ) 



belongs to A^ +2 . By calculus rules in the Fourier space, since F is smooth, it 
holds 

(F{F + h 2 F^) 



vM(x,v)dv = iv 



[x, 0). 



v \ v — i5V(x, rf) 
By performing the derivative and then taking into account that 

T{F + fr 2 F (2) )(x, 0) = FF(0) = 1, V„.F(F + fr 2 F (2) )(x, 0) = 



and that (V v 6V)(x,0) = V x V(x)/m, one gets (Oil . 

Analogously, the second moments of M are well defined, and, by calculus rules, 
it holds 



ViVjM(x, v) dv 



d 2 f(F + h 2 FW) 
drjidrjj v — i5V 



(x,0), \/i,j = 1,. . .d. 



11 



and 



d 2 f TF \ 1 (&TF\ Q] 1 / ggV ggVA 

drjidrjj \v — iSV J ' z/ \drjidrij J ' z^ 3 \ <9r/j c^fo / 

1 2 <9V(x) 



z//3m z/ 3 m 2 c?Xj dx 



J 



d 2 ( tfw \ = i f ^re^ N = g gM^) 

Thus the thesis follows. ■ 

Remark 4.2 The state that describes the system under the effect of the interac- 
tion with the environment and of the strong potential is described by the function 
nM, with M defined by (14. 2p . The fluid velocity relative to such state is non- 
null and given by (14. 8p . In contrast, the velocity of the system in the state w eq 
defined by (12. 7\\ (i.e., when it is subject to the sole influence of the environment), 
is J vw eq dv = 0, as expected since it is an equilibrium state. Moreover, the 
expression of the second moment tensor (14. 9ft has to be compared with 

v <g) v w cq dv = n[ h — — -V (g> W 

\pm 12m z 



They differ by the second summand in (14.91) that is to be referred to the strong- 



field assumption (cf. [32]). □ 



As a consequence of Proposition 14.11 the following subspace is well-defined 

(X k ) M := {a(x)M(x,v),a E L 2 X } C X k , 

which coincides with ker(^4 + C) when A + C is considered as an operator on X k . 
Accordingly, we can decompose the space X k as 

X k = (X k ) M © (X k )° (4.10) 

with 

w(x, v) dv = 



(X k ) u := E X k 

and define the corresponding spectral projection V from X k into (X k ) M , by 

Vw := M / w(x,v)dv, 



M 



and Q := X — V . The following corollary is still a preliminary result for our 
asymptotic procedure. 
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Corollary 4.1 Let V G with d- admissible k. Then, the operator Q(A + C)Q 
is an isomorphism of [Xk)° onto itself, with 

\\A + C\\ B{Xk) <C{d,k){\\V\\ m + v). (4.11) 

If, in addition, V G H* +j with j > 0, then Q(A + C)Q is an isomorphism of 
(Hi) onto itself, with 

\\A + C\\ B(Hi) < C(d,k,j) (\\v\\ H ^ +v) • ( 4 - 12 ) 

Proof. The operator Q(A + C)Q, when considered as an operator acting on 
(X k )°, reduces to 

Q(A + C)Qu = Q[V]u-vu, \/ue(X k )°. (4.13) 

Then the thesis follows from Lemma |3"TT| Prop. 14.11 and by the skew-simmetry of 
the pseudo-differential operator. The second statement and estimate (14.121) can 
be proved analogously. ■ 



5 The asymptotic expansion 

According to the decomposition of the space Xk, every function w G X k can be 
written as w = Vw + Qw, with Vw G (Xkju and Qw G (Xk)° . Let us call 
if := Vw and if) := Qw. Observe that, for all w G X kl f Vw(x,v) dv = n[w](x), 
while f Qw(x, v) dv = 0, that is, we separate the part of w that contributes to the 
density n[w] from the other one. Precisely, it holds Vw = n[w]M, by definition. 

Applying formally the projection V, respectively Q, to the Wigner-BGK equa- 
tion (13.51) with unknown w, we obtain the following system of equations with 
unknown <p and if) 

= VSV^ + VSQif) 

1 (5-1) 
= QSV<p+ QSQif)+ -Q{A + C)Qif) 

where we used (A + C)V(p = and V(A + C)Qip = 0, together with the initial 
conditions 

p(0) = = Vw , ^(0) = = Qw . (5.2) 

System (15.11) consists of an evolution problem with unknown functions tp = n[w]M 
and if), and it is supplemented by the initial conditions (15. 2ft . It is a reformulation 
of (1331). 

Since we expect the solution w to be subject to rapid changes for small times, we 
split the functions ip and if) into the sums of the "bulk" parts ip and if) and of the 
"initial layer" parts <p and if), 

<p(t) = (p{t) + , = ip(t) + ipr- 
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The bulk part <p is left unexpanded and the other parts are expanded in terms of 
e as follows 



(p{r) 

m 



(p {r) + e(pi{r) + e 2 (p 2 {r) + . . . 

+ e$i(t) + e 2 4) 2 {i) + . . . 
Mr) + e^i(r) + e 2 ^ 2 (r) + . . . 



(5.3) 



with t = t/e. Accordingly, Eqs. (15.11) for the bulk part terms of the expansion up 
to the order e 2 become 




VSVip + VSQtPo + eVSQ^ 

Q(A + C)Q^ 

QSV<p + Q(A + C)Q$ X 



(5.4) 



while the equations for the initial layer parts read 
( dfio 



Or 
dip i 

dr 
dr 
dr 



0. 

VSQMt) 

Q(A + C)QMr) 
Q(A + C)QMt) 



(5.5) 



and the initial conditions ( 15. 2ft yield 



^(0) + ^ (0) + e^(0) 
^o(0) + ^ (0) 
&(0)+&(0) 



QSQMr) 



Po 
0. 



(5.6) 



System (15.41) . together with fj5.5H - fj5.6p . is an (9(e 2 )-approximated version of (15.11) 
with (15. 2p . once the expansion (15. 3p has been introduced. In fact, the equations 
in ( 15. 4p can be decoupled: by Corollary 14.11 the operator Q(A + C)Q is invertible 
in (X k )°, thus 







{Q(A + C)Q)' 1 QSV<p 



which implies 



dip 
~dt 



VSVip - eVSQ(Q(A + C)Q)- l QSVp . 



(5.7) 
(5.8) 



(5.9) 



Thus, system (15.41) reduces to the system fl5.8D - fl5.9p . with unknown functions 
(p(x, v, t) = n(x, t) M(x, v) and ^)\. Next section shall be dedicated to reformulate 



14 



Eq. (15.91) as an equation with unknown n. The analysis of system (15. 5p . with 
unknown if and ip and initial conditions (15.61) . is postponed to Section [3 it shall 
provide an appropriate initial condition for Eq. (15. 9p . Finally, in Sections 7 and 
8 we shall establish a well-posedness result for the approximated problem. In our 
main theorem (cf. Thm. 19.11) . we shall prove that the solution ip + ip of equations 
flE0 indeed differs from [(p(t) + <p (r) + e£i(r)] + $ (*) + <$i(*) +$o{r) + e$i(r)], 
satisfying the approximated problem (I5.4p - (I5.6I) . by a term of order e 2 . 



6 The high-field quantum drift-diffusion equa- 
tion 

The aim of the present section is the reformulation of the abstract equation (15.91) 
as an equation with unknown n. 

Lemma 6.1 Let V G H^ +2 with d- admissible k . Eq. (15.91) with unknown tp(x,v,t) - 
n(x,t)M(x,v) can be rewritten as an evolution equation with unknown n(x,t) of 
the form 

- V- (nW) V- Vn 

at vm upm 

[v- (n(v ® v)y vy) + v- v- (nvy ® vy)] 



is 3 m 2 



eBft 2 

~ -7PrV-V-(nV(8)Vy) = (6.1) 

Remark 6.1 The first line of (16. ip consists of the terms of the classical DD 
equation. The second line is peculiar of the strong-field assumption, being a 
correction of order e, and consists of the additional term 

1 W (8) VV 
v v 2 m 2 

in the pressure tensor, and of the term 

1 (V®W)VWV\ 



f v®v)vvv \ 

\ i> 2 m? J 



in the drift term. Both terms are quadratic in the potential V. The second line 
can also be written as 

V- [W <g> VVVn + n (2V ® WW + AVW)} . 



is 3 m 2 



This expression is the same obtained in [32] from the semi-classical Boltzmann 
equation with high-field scaling. The last line is the quantum pressure term 

(cf. usiizil). □ 

The proof requires the following preliminary lemmata. 
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Lemma 6.2 Let V £ H^ +2 with d- admissible k. Then the equation 

(A + C)w = M t-v+ vMdv) , (6.2) 



admits a unique solution (D 2 )i £ (Xk+i)° , Vi = 1 , . . .d. Moreover, let D be the 
matrix defined by 

D i:j (x) := J v i (D 2 ) j (x,v)dv, 

then 2 

D ( x ) = - \7T- + ® W + ^T^ V ® W l ( x ) • ( 6 - 3 ) 

z/ \pm z/ 2 rrr 12m^ / 

Proof. Since the right-hand side of Eq. (16.21) belongs to (Xk+i)°, it satisfies 
the compatibility condition (14.31) and there exist (D 2 )i £ (X^+i) , Vz = 1 , . . .d 
satisfying ( 16. 2ft . More explicitly, D 2 solves 



{Q[V} - v)D 2 (x,v) = -M{x,v) 



VV(x) 
v + ^ 



z/m 



(6.4) 



since j D 2 (v)dv = and by (14.8| . Multiplying the left-hand side of Eq. ( 16.41) by 
v + VVy (z/m) and integrating over IR d we obtain 

v+ W ( x ) ) g) [e[y] _ „] £) 2 ( X) = _^ [ Vi (D 2 )j(x,v)dv, 
z/m / y 

by using the skew-simmetry and (D 2 )i £ (X fc+1 )°. Thus Eq. (16.41) gives 

W(x)\ / W(x) 



uDM = (v + -^^)®(v+^^)M{x,v)dv (6.5) 
J \ z/m / V vm J 

v ® uM(z, v) cfc — Vyfs) <g> Wfz) . 

z/^m 2 



From (14. 9p . the thesis follows. ■ 

Remark 6.2 By considering the expression (14. 8p for the fluid velocity, we can 
recognize in (16. 5p the classical definition of the pressure tensor in terms of M. This 
is to be expected, since the function M is the solution of the the evolution problem 
with e = 0. Thus, the term with diffusion tensor D is what we expected to find as 
correction of first order in e. By (16.31) . it consists of the standard temperature and 
quantum pressure tensors, and of the additional tensor l/(u 3 m 2 )'VV(x)<E)'VV(x), 
to be referred to the strong-field assumption (cf. Remark 14. 2p . □ 

Lemma 6.3 Let V £ H^ +2 with d-admissible k. The following equation 

(A + C)w = -v- W X M + M Jv W X M dv , (6.6) 
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admits a unique solution D\ £ (Xk + i)° . Moreover, the vector W defined by 

W(x) := J vDi(x,v) dv 

can be calculated explicitly 

W W = - ( 2 Zf^W (x) + ^ + ^V- V 8 w) W . (6.7) 

Proof. Under the regularity assumptions on V, M £ -ffjj+i := H\ ® and 
the right-hand side of Eq. ( 16. 6ft belongs to (X^+i) , thus it exists -Di £ (X^+i) 
solving 

(6[V] - u) = -v V x M(x, v) - M(x, v)V- ^ 



which is equivalent to (16. 6p . since J Di(v) dv = 0. Multiplying by v + W/um 
and integrating over H d we obtain 

v JvD 1 (x,v)dv = J(v+^^jV x -(v + ^^jM(x,v)dv 

= Vj.- \v®vM(x,v)dv — — AWx)W(x). (6.8) 

The thesis follows directly from (14 .9p and 

v- (w ® w) = avw + (v ® v) v w . 

■ 

Proof of Lemma 6.1 First of all, let us write the explicit expression of the 
operators appearing in Eq. (15. 9p . Observe that V £ H^ +2 implies M £ H\ +l . By 
definition, 

VSV(p = -M JvV x (nM)dv = - MV X - (^n jvMdv^j , 

and (VSQg) = (VSg) - (VSVg), i.e., explicitly, 

(VSQg) = -M (JvV x gdv-Jv x gdv-JvMdv-JgdvV x -JvMdv\ . 

(6.9) 

Moreover QSVip = (SV — VSV)ip, i.e., explicitly, 



QSVCp = n 



-vV x M + M v-V x Mdv 



+ V x n- 



M [ —v + 



JvMdv^j 



. (6.10) 



By Lemmata 16.21 and 16.31 D 2 (x,v) = (D 2 )i(x,v) and Di(x,v) are solutions with 
(D 2 )i,D 1 £ (X k+1 )° of Eqs. (16. 2p and (16.61) . respectively. Then, by some manip- 
ulations, 

VSQ(Q(A + C)Q)- 1 QSV<p = VS{D T Vn) + VS{D 1 n) , (6.11) 
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where the right-hand side can be written explicitly as 



VS(D 2 -Vn) = —M J v- V X (D 2 - Vn) dv = -MV 3 



v (g> D 2 dv ) • Vn 



VS{D in ) 



-M J v- V x (Din) dv 



■MVx- 



n I vD\ dv 



Hence, by simplifying the common factor M, Eq. (15.91) reads 

dn 
~dt 



- V*- (^n JvMdv^j + eV x - (D- Vn + nW) , 



and the thesis follows by using f!6.3p and (16.71) . ■ 
As a consequence of (16.111) . Eq. (15.81) defining the other non-zero term of the bulk 
part expansion, ifii, can be rewritten as 



V>iO,t) = -RVn + D in } . 
The explicit version of this expression shall be given in Eq. (18.201) . 



(6.12) 



7 Rigorous results: the initial layer part 

The aim of the present section is to prove existence and regularity of the solutions 
of Eqs. (15.51) . together with the initial conditions (15.61) . The first equation in (15.51) 
yields 

0o(r) = 0, 

since we expect that lim^oo (po{r) = 0. The equation for tpo with the appropriate 
initial condition coming from (I5.6p - (I5.7I) . is 



dipo 
dr 
MO) 



Q(A + C)Qi; 



(7.1) 



We recall that the operator Q(A + C)Q on (Xk)° reduces to 
Q(A + C)Qw = Q[V}w-vw 1 Vw6(IJ 



(cf. (14.131) ). By the product shape in Fourier- variables of the pseudo-differential 
operator (cf. (12. 3p ). it is more convenient to consider the equation for Fipo, that 
looks like 

d ~ ~ 
— Fil>o(x, r], t) = (i 5V(x, 7]) - vjFipoix, r], t) . 

Thus, we define, for all w G L 2 (H 6 ; 1R), the semigroup G(r) 

G{t)w(x,v) := T~ x (e {iSV{x ^- u)T ^w(x,r])) (7.2) 
= e~ VT T~ l {e iSv ^ T Tw(x, 77)) , Vr > . 

The function ipo(r) = G{r)ipQ formally satisfies system (17.11) . Moreover, 
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Lemma 7.1 If w G (X^) and V G iif* with d- admissible k, then there exist 
< Uk < 1, and a constant C(\\V\\ H k) > 0, such that 

\\G(T)w\\ Xk <C(\\V\\ m )e~^\\w\\ Xh . (7.3) 
If, in addition, w G {H^f an d V H^. +: > , then 

\\G(t)w\\ h{ < C{\\V\\ Hi+j ) e-^\\w\\ Hi , (7.4) 

with appropriate C(\\V\\ H k+j) > 0, and < v>k+j < 1- -So- (17.21) defines a strongly 
continuous semigroup on {X^) (respectively on (Hi) ). 

Proof. By definition we have 

\\G(r)w\\ Xk < Ce-^(\\e i6V ^Fw(x, V )\\ L ^ + \\V k v (e iSV ^Fw(x, V )\\ Li ^ 

< Ce— {\\w\\ Llv +P k (r\\V\\ m )\\w\\ x ^ 

< e-^max{e-^ T P k (r\\V\\ m )}\\w\\ Xk , 

r>0 

where < v k < v and P\~ is a polynomial of degree k. The estimate (I7.4p can 
be proved analogously. The last assertion follows immediately by applying Hille- 
Yosida Thm., thanks to (I7.3P (respectively (I7.4p ). ■ 

With Lemmata 14.11 and 17.11 we can prove the following proposition. 

Proposition 7.1 // w G H\ +l and V G H^ +2 , with d-admissible k, then all 
terms of the initial layer expansion are well-defined and satisfy the following es- 
timates: 

\\Mr)\\ Xk < M ie -^|Klk, (7.5) 
\\0i(r)\\ Xk < M 2 e-^\\w \\ Hl+i , (7.6) 

\\Mr)\\ Xk < M 3 e-^|KIUi +1 , (7.7) 
for some constants Mi, M2 and M3 (depending on norms ofV). 
Proof. The unique solution of system (17. ip is 

Mr) = G(r)^ , (7.8) 

and (17. 5p follows immediately from (17. 3p since ipo = Vwq G X^. Now we shall 
consider, among Eqs. (15.51) . the following one: 



Ot 



-(t)=VSQMt). 
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The right hand side is well-defined by considering the definition of the operator 
VSQ (cf. (JS3D), together with Lemma EH since ^ G and V G #£' +2 . 

By integrating with respect to r and considering lim r _,. 00 (p\{r) = 0, we obtain 

/oo 
VSQfa{s) rfs 

P5Q[Q(^l + C)Q]- 1 [Q{A + C)Q]G(s)^o ^ = 

/oo 
Q(^ + C)QC?(3)^o da- 
Last integral is well-defined since the integrand Q(A + C)QG(s)fa is equal to 
G(s)Q(A + C)Qfa, which is continuous in H\ +l (by Lemma I7TTI) . Moreover VSQ 
[QiA + tyQ]- 1 G C(Hl +1 ,X k ) (by Corollary gi]), then it can be taken outside 
the integral. Since 

thanks to the exponential decay of G in H\ +l and the continuity of the operator 
VSQ[Q(A + C)Q}~\ we obtain 

0i(r) =VSQ[Q(A + C)Q]~ 1 G(r)ij , (7.9) 

and, in particular, 

<fi(0)=VSQ[Q(A + C)Q]- 1 fa, 
which provides the initial datum. Then, (17.61) follows from the estimate 

Wi(r)\\x k < \\\VSQ[Q{A + C)Q]- l \\\ ||G(r)^||^ +1 < ( 7 - 10 ) 
< M 2 e-^ T \\M\m 



fc+i 



where ||| • ||| denotes the norm in C(H\ +l , Finally, we prove that the equation 

^l( r ) = Q(A + C)Qfa(r) + QSQfa{r) 
or 

is classically solvable. The initial condition for ipi can be obtained from Eqs. 05. 6|) . 
together with Eq. (15.81) for fa, 

fa{0) = -fa(0) = [Q(A + C)Q]- l QSVm, 

and by considering 

£(0) = <A) - £o(0) - e0i (0) = - eVSQ[Q(A + C)Q]~ l fa. (7.11) 

Since ^i(O) is by itself a correction of order e, we neglect the term of order e in 
the expression for <^(0), and it yields 



^i(O) = [Q{A + C)Q]- l QSV<p 



o • 
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By Lemma I7TTI G is a semigroup on (Hl +1 )° and, thanks to the assumption on wq, 
ipo is in the domain of Q(A + C) Q when it is defined in D(S) D (Hl +1 )°. Therefore 
V , o( r ) = G(r)ipo is differentiable on [0, oof in X^+i so that the inhomogeneous 
term QSQipo(r) is differentiable on [0, <x>[ in X^. This, together with V'i(O) = 
(QiA + tyQ^QSPtpo e D{Q{A + C)Q), shows that 

Mt) = G(t)M0)+ f G(r-a)QSQG(a)ij da (7.12) 

is a classical solution. The estimate (17. 7p follows from [Q(*4 + C)Q] -1 QtS"P G 
£(Hl,X k ) and from (17^ : 

< K ie -^||^o|U fc +K 2 e-^ /%^-^||^o|ki^ 

./o 

< K 1 e-^ T || ¥ »o||x fc + K 3 e-^+ aT ||V'o||^ 1 < M 3 e- v ^ T \\w \\ H i . 



In order to obtain an initial value for Eq. (16. ip with unknown </? = nM, we 
consider again (17.111) . Let us call n (x) = f w (x, v)dv, such that (p = n M 
and, by dividing both sides of the expression (17.111) by M, it yields 

n(x,0) =n (x) + e Jv V^ 1 ) fr, v) dv , (7.13) 

by using the explicit expression the operator (Q(A + C)Q)~ l (cf. ( 14.131) ). In the 
following we shall call 

n(x, 0) = n (x) + enAx) with nUx) := [v-V x J-'~ 1 ( — — - — J (x, u) dv . 

J \idV-vJ 

(7.14) 

The explicit expression for (17. 9p can be obtained analogously and reads 

8 Well-posedness of the high-field QDD equa- 
tion 

In this section, we establish a well-posedness and regularity result for Eq. (16. ip . 
with a given external potential V. The equation can be rewritten in divergence 
form as 

din 

— -Vn-Gn-£n = 0, (8.1) 
dt v ' 

where we indicate 

Vn = eV ■ (DVn) , Qn = eV ■ (\N n) , 5n = V-(En) 
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with 



D = D(x) := - — + — -W®W + -^V®W)(x) 

W.W W := I( 2 X|ZI W+ ^ + -f_V.V»V^)(x) 

v \ u z m z u z m z 12m z 

E = E(z) := V*GlR rf . 



Assumption 8.1 V belongs to H x +2 with d- admissible k and it satisfies the fol- 
lowing 

3c>0 s.t. D(x)y®y > c\y\ 2 , \fx,yeJR d , 

This implies that D is a uniformly elliptic differential operator. Thus we can 
state the following: 



-l,oo 

X 

2 



Proposition 8.1 LetV satisfy Assumption \8. II and, in addition, VAV^ G W 3 
with j G IN. Then the unique global solution n = n(t) of Eq. (18.11) with n(0) G L 
satisfies n(t) G H 3 X for t > 0, and the following estimate 

\\n(t)\\ Hi <M,(et)^/ 2 \\n(0)\\Ll ( 8 - 2 ) 
holds with Mj > 0, for e, t — » + . 

In the following, by VF G L 2 (IRf) we mean VF G (L 2 ) d . Moreover, we consider 
< e < 1 and the constants are independent of e, unless specified. 

Proof. By Assumption l8.1l on the potential V, the operator D defined on D(T>) = 
H 2 generates an analytic contraction semigroup (T(t))t>o on L 2 (cf. Thm.VI.5.22 

of p3]). 

Let us derive here a basic estimate we shall use intensively in the following. By 
Assumption 18. 1[ for all u G H 2 

|| V- Vllq < C||DV « Vu hi < § ( Mis + ||divD. VHU.) , 

where the second inequality follows from the definition of the operator T>. More- 
over, by using that for all u G H 2 , 

\\Vu\\ L 2 <C s \\u\\ L 2 +5\\V-Vu\\ Li , W5>0, (8.3) 

with C s > 0, it holds 

||divD- Vu\\ Ll < ||divD|| ig =(as||M||^ + 5\\V- Vu\\ L 2) , V5 > 0, 
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with C$ > 0. Thus, 

l|V ' V " h i £ t (l^-||d.vD l U.A) (l|dlYD|U - Cill " IU ' + ' 

and, in conclusion, for an appropriate choice of 5 > 0, exists some constant C > 
such that 

HV-Vullig < -(NUg + H^lUg). (8.4) 

The operator Q can be written as Q = Q\ + Q2 with ^1/ := eW- V/ defined on 
if*, and £? 2 / : = £ V- W/, defined on L 2 .. The operator £1 is P-bounded, i.e., for 
all n G £>(£>), 

HMUl < e||W|| iso C(||n|| Li + ||V-Vn|| Li ) (8.5) 

< ||W|| iS0 (C||n|| Li +C||X>n|| Li ) 

< 6||n|| L | + a||Pn|| L | , (8.6) 

by using (18. 4p and e < 1. Moreover, the D-bound a x defined by 

ai := inf {a > | 3 6 > s.t. ([HID holds} 

is zero, by substituting (I8.5P with fl8.3j) . The operator C/ 2 / : = eV- W/, defined on 
L x , is bounded. The operator £ can be written as £ = Si + £ 2 with E\f := i?- V/, 
defined on if*, D-bounded with P-bound a 2 = 0, since, for all n G D(T>), 

< C||£|U« U7f||ra|U 2 + ^(N^ + IIPnll^V V5>0,(8.7) 

by using (18.31) and (18.41) . The operator £ 2 / := V-i?/ is defined on L 2 and 
bounded. 

Thus, by Thm. III. 2. 10 of [16], (£> + (? + £, D(T>)) generates an analytic semigroup 
on L 2 X that we shall indicate with (S(t)) t > . More precisely, it holds 

\\(V + g + S) a S{t)u\\ Ll < M a r a \\ U \\ Ll , t -> + , Va > (8.8) 

with M Q independent of e, by employing Lemma III. 2. 6 of [T6] . 

In order to derive estimate (18.21) . let us start from the following inequality 

e m/2 \\u\\ H ™ <C\\(V + g + £) m ' 2 u\\ Ll (8.9) 

that is yielded by similar arguments to (]8.4j) . By combining (I8.9[) with (18.81) . we 

get 

llS^MIiri <Ce-^ 2 ||(P + ^ + 5)^(^11^ <C,(et)^/2|| n || Li; (8 . 10 ) 

which holds for all u G L\ and for small t. ■ 
In estimate (18.21) can be easily removed the singularity with respect to t: 
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Corollary 8.1 Let the assumptions of the Prop. (18. ip hold. In addition, let n(0) 
belong to H x . Then the solution n(t) belongs to H x for allt>0 and satisfies 

\\n(t)\\ Hi <Ce-^\\n(0)\\ H i, (8.11) 

for e, t^ 0+ . 
Proof. 

The following inequality holds 

e m/2 |Mk™ < C\\(V + g + £) m / 2 u\\ L 2 < C\\u\\ H rn (8.12) 

for all u G D((V + Q + £) m/2 ), with m < j, and for e 0+ (cf. Then, in 

particular, 

ii^Hh™ <c^ m /ioD+£+£r/ 2 s(t) M |i^ 

< Ce- m / 2 ||(I? + ^ + £) m / 2 M || L 2 <Ce- m / 2 || M ||^, (8.13) 

which holds for all u G D((V + Q + £) m / 2 ), for small t and e: the first inequality 
sign corresponds to the first one in Eq. (18.121) . the second inequality follows by 
exchanging (T> + Q + £) m / 2 with S(t) and the third one comes from Eq. (I8.12p . ■ 



Remark 8.1 Observe that in the low- field case, the QDD equation looks like 
Eq. (JHIIJ with 

W(s) = ( t|^V- V ® VV) (x) , D(x) = -(^- + ® W ) (*) • 

In order to establish well-posedness result and estimate (18.21) for all j G IN, the 
same assumptions of Prop. 18.11 are required on the potential V and on the initial 
datum n(0). This result is to be compared with the analysis in [21], where it 
is tackled the fourth-order, non-linear equation obtained by the approximation 
Vlogn = -(3VV + 0(H 2 ), cf. [19]. D 

By increasing the assumptions on the initial datum, we can remove the singular 
behaviour of the estimate (I8.2p with respect to t and e. 

Corollary 8.2 Let V satisfy Assumption {8J\ and VAV G W 2j ~ h °°. Then the 
solution n(t) of Eq. (18.11) with n(0) G D{V 3 ) satisfies, for e, t — > + 

\\n(t)\\ Hi <C\\n(0)\\ H *j. (8.14) 

Moreover, the following refinement holds 

H*)ll^<C||n(0)||^. (8.15) 
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Proof. We prove the thesis in the case j — 1. For j > 1 the thesis follows by 
an induction procedure similarly to [8]. Due to the regularity with respect to the 
variable x of the solution n(t), for t > 0, we can find the evolution equation for 
Vn by differentiating 

V \^nj = ^(Vn) = V(P + Q + £)n = (P + Q + 5) Vn -[(V + g + E), V]n, 

where we indicate with [(V + g+E), V] the commutator among the two operators. 
Since 

-[(V+g+E), d k ] = eJ2d* 10,D, ,<),„) ■ ^\y, (d k (eW ! + E*) n) =: (V'+g'+E')n, 

hi i 

Vn satisfies 

A(Vn) = (X? + £ + £)Vn + (!?' + £' + £')«• (8.16) 

The solution of the previous equation can be expressed by the Duhamel formula 
via the analytic semigroup S(t) generated by (V + Q + E), as 

Vn(t) = S(t)Vn(0) + [ S(t- s)(V + Q' + E')n(s) ds . (8.17) 

Jo 

Moreover we can estimate 

\\Vn(t)\\ Li < C\\Vn(0)\\ Ll + C [\\(V' + g' + E')n(s)\\ Ll ds 

Jo 

< C\\Vn(0)\\ Ll + C f(e\\n(s)\\ m + \\n(s)\\ H i)ds (8.18) 



< C\\n(0)\\ H , + C / \\Vn(s)\\ L2 Js. 

Jo 

by using ( 18. 1 If) with j = 2, provided n(0) G Therefore, by Gronwall lemma, 
we derive (18.14p . In order to prove ( 18.151) . we apply for the function n(t) = 
S(t)n(0) the first inequality in (18.101) with j = 2 and we obtain 

\\n(t)\\ H 2<^\\(V + g + E)S(t)n(0)\\ Li . 

Then we use (E2D for the term (V + Q + £) l / 2 S(t) [(V + Q + £) 1/2 ™(0)] with 
j = 1, and we get 

\Ht)\\m < ^Wiv + g + Ef^s^v + g + Ef'MmLi 

< -t~^\\{v + g + EfiMmLi 

e 

< £ r l/2|| n(0) || (8 . 19) 

e 
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where for the last inequality the estimate (18.121) with m = 1 is used. Hence it 
holds for all n(0) G V((V + Q + Efl 2 ) = E\. By using (JHHSD in (JHUBJ, we get 
( 18.151) in the case j = 1. ■ 

The other (non-null) term of the bulk part is ipi, which is of first order in e. Since 
it satisfies 

$i = -{Q(A + C)QY l QSV(nM) 

(cf. Eq. (15.81) ). by using the definitions (16.101) and (14.131) . it can be written ex- 
plicitly as 

v, t) = Vn(x, t)- T { — T vM + } (x, v)+ 

{to V — v \ vm J J 

n(x, t) T~ x { — )- T (v- W X M + — - ] 1 (x, v) , V fx, v, t) G 1R M x 1R+ . 

\io V — v \ vm J J 

(8.20) 

Thus, the estimates for the solution n in the previous corollary are the crucial 
ingredient to establish well-posedness of the definition (I8.20p and the behaviour 
with respect to time of the function and of its derivatives. 
Another fundamental aspect is the shape of the initial datum n(0) for Eq. (18. ip : 
by (I7.14p . it is given by n(0) = n + erii, and the following estimate holds 

||n(0)|U < ||no|U+e||ni||^ < \\w \\ H , + eC(\\V\\ H k + , + 2)\\w \\ H j + i, (8.21) 

^ k-\-l 



for all ^-admissible k, by using the estimate (I4.12p . 

Proposition 8.2 Let n be a solution of the drift- diffusion ( 18.11) with initial value 
n(0) given by (17.141) . with wq G H^ +1 , and with V satisfying Assumption \8.1\ and 
VAV G W^'°° . Then ip\ is strongly differentiate with respect to t > 0, and for 
every t > it satisfies 

Mt) eV(Q(A + C)Q)nV(QSQ). 

Moreover there exists a constant M > such that, for e, t — > + , 

(8.22) 
(8.23) 
(8.24) 

Proof. If we differentiate with respect to t the expression (18.201) . the only t- 
dependent functions are Vn and n, explicitly 

d t ipi(x,v,t) = d t (Vn)(x,t)- A(x, v) + d t n(x, t) B(x, v) , V(x,v,t) G R M x 1R+ 

(8.25) 



ptMt)\ 




< 


M \\wo\ 


k 4 +i ' 


\\dtMt)\ 


k 1 


< 


M(l + 


l/t)\\w \\ H 4 +i 


\\SQMt)\ 




< 


M||w | 


1^+1 • 
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where the functions Ai,B, defined in (18.201) . are sufficiently regular, because of 
the assumptions on V. The differentiability of n with respect to t depends on 
the analiticity of the semigroup (S(t)) t >o- The differentiability of Vn, instead, 
follows from the expression (18.1 7ft : since each term is continuously differentiable 
in time, also Vn is. 

Moreover, by using <9 t Vn = Vd t n and the evolution equation for n, 

\\d t (Vn)\\ Li = \\V(V + g + S)S(t)n(0)\\ L 2 

< \\(V + g + £)S(t)n \\ H i + e\\(V + Q + £)S(t)m\\ H i 

< \\S(t)n \\ H 3 + e\\S(t)ni\\ H 3 < C(\\w \\ H s + e\\w \\ H i +i ) 

< C\\w \\h Ui (8.26) 

where we split n(0) = no + eni and we use ( I8.15p . together with (18.211) . Similarly, 
\\d t (Vn) \\ Hi < C (|| V(V + g + S)S(t)n{0) \\ Li + || V- V(V + g + S)S(t)n(0) \\ Ll ) 
< C [\\w Q \\ Ht+i + \\{V + g + £)n \\ m + L\\ {V + g + £) ni \\ L ^ 



< C^\\wo\\ Ht+i + \\wo\\ Ht + -\\w \\ Hl+i ) , (8.27) 

where the first addendum in the inequality comes from (18.261) . The second and the 
third terms come by exchanging S(t) with (D + g + £) and using n(0) = n + erii, 
then we apply estimate ( 18 . 1 5f) to get the second term, and estimate ( 18. 2p to obtain 
the third term. Finally, inequality ( 18.271) follows from estimates ( 18. 12ft and ( I8.2ip . 
In order to prove (I8.24[) . let us consider again the abstract definition of z/>i (see 

<my. 

Mt) = -(Q(A + C)QY\QSV)(n(t)M) . 
Since QSV(p reads (see ( 16.101) ) 



QSVnM = n 



-v-V x M + M v-V x Mdv 



V x n- 



M [ —v + vMdv 



under the present hypotheses, QSV{nM) belongs to H%, thus ipi(t) e T>(Q(A + 
C)Q)nV(QSQ). By (ESQ}, 

SQipi(x, v, t) = V- Vn(i, t) v ■ A(x, v) 

+ Vn(x, t)-(v-V-A + vB)(x, v) + n(x, t) v ■ VB(x, v) . 

Thus, in order to estimate HSQ^i^Hxfc and ||<S'Q'0i(t)||_H-i, it is necessary to eval- 
uate || V- Vn(t)|| £ 2 and ||V- Vn(t)\\ H i, || Vn(t)|| L 2 and ||Vn(t)||#i, respectively. In 
particular, 



|V- Vn{t)\\m < C\\S(t)n(0)\\m < C7||ra(0) < C\\w \\ H s + eC\\w \\ Ht 



fc+i 



again by (I8.15p . Thus, we can conclude by using the regularity properties of 
A U B. m 
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Remark 8.2 Observe that it is possible to remove the singularity for t — > + 
in the estimate (18.231) . by assuming wq G H^ +1 and modifying last two lines of 
(E2ZD as follows 

\\d t (Vn)\\ H i < C (\\w Q \\ Ht ^ + \\{V + g + £)n \\ m + e|| (£> + £ + ^mH^) 
< C (lko||j?4 +1 + ||wolU* + e||«;o||H« +1 ) ■ ( 8 - 28 ) 

□ 



9 Estimate of the error 



In this section we prove rigorously that the high-field QDD equation, originated 
by the asymptotic expansion up to the first order in e, is an approximation of 
order e 2 of the high- field Wigner-BGK system (13. 51) . To this aim, we consider the 
errors obtained by replacing the functions Vw = ip and Qw = ip by the terms of 
their expansion up to first order in e. We shall prove the following 

Theorem 9.1 If the initial value wo belongs to H^ +2 and V satisfies Assumption 
lg.il and VAV" G W^°° , then for any T, < T < oo, there is a constant C 
independent of e such that 

<p(t) + ip(t) - [<p(t) + e^i(t) + Mt/e) + eCp^t/e) + e^it/e)] < Ce 2 , (9.1) 



uniformly for < t < T. 



This result relies on the estimates established in Propositions 17.11 18.21 about the 
behaviour with respect to time of the initial layer functions (f\, ip\ and the bulk 
functions. Let us split the error in two contributions 

y(t) = ipit) - [(p{t) + e£i(r)] , z(t) = - $ (r) + e^(t) + e^(r)] (9.2) 

where r = -. The evolution equations for the errors y and z can be deduced 
from those satisfied by their components (cf. systems (15. II) . (15. 41) . (15.51) ). Hence, 
we have 

% = VSVy + VSQz + f 

Z i (9 - 3) 

_ = QSVy+ QSQz + -Q(A + C)Qz + g 



dt 

with initial conditions 



0, 



2/(0) =0, z(0) 
and inhomogeneous terms / and g defined by 

f(t) = e [vSP<P!(t)+VSQMt) 

git) = e -^ + QSQMt) + QSV^ 1 (r) + QSQMr) 
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It is convenient to separate the evolution of the error relative to initial layer part 
from the one corresponding to the bulk part. Let us define 

r = y + z = ri + r b 

with 

Ti = -e0! - i) Q - eV>i , r b = y? + ip - - erfi . (9.4) 

The derivation of estimate (19.11) is split according to (19.41) in the next two Lem- 
mata. 

Lemma 9.1 Under the assumptions V G H k+A and G H% +2 , for any T, 
< T < oo , there is a constant C independent of e such that 

\\ri{t)\\ Xk <Ce\ (9.5) 

uniformly for < t < T . 

Proof. The initial layer error n = ri(t) satisfies the equation 
Br 1 ~ ft\ 

~^{t) = Sn(t) + -(A + C)n{t) + eS(0! + Vi) I- J , n(0) = . (9.6) 

The operator S + (A + C)/e generates an uniformly bounded semigroup in X k , 
Z(t), cf. [IB]. Thus, the mild solution of (19.61) is given by 

nit) = Z{t)n(0) +e [ Z{t- 8)3(0! + 4>!)(s/e)ds , 







with 



\\ri(t)\\ Xk < Ce [ \\S(0! + i h )(s/e)\\ Xk ds. (9.7) 
Jo 

The estimate of HS^i+^i) (s/e) \\ Xk is a bit tedious, thus we simply sketch it. It 
is convenient to use the projections V, Q and evaluate VSV0!, QSV0!, VSQipi, 
and QSQip! separately. By their definitions (cf. (16T9]) . (16TT01) ) VSV, QSV G 
C(Hl,Xk), provided V G H^ +2 , thus, it holds the following modification of the 
estimate (I7.10p for X , 

\\VSV0!(r)\\ Xk < IWVSVWWWVSQIQiA + VQr'lW ||C(r)^ ||^ +1 

< Me-^ T ||Vo||^ +1 , (9-8) 

since VSQlQiA + QQ}" 1 G C(H 2 k+1) Hl), provided V G H k x +i . Analogously, 

\\QSV0!{r)\\ Xk < \\\QSV\\\\\\VSQ[Q(A + C)Q]- l \\\ \\G(r)M\ Hl+i 

< Me-^HVoll^. (9.9) 
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Let us recall the expression for ipi (cf. ( I7.12p ) 

Mr) = G(t)&(0) + / G(t- a)QSQG(a)i> da 



We evaluate WVSQ^t)^ and WQSQ^t)^. Both PS 1 Q and QSQ belong to 
C(Hl +1 ,X k ). Moreover [Q(A + C^^QSP G C(H\,X k ) by definition, provided 
V G H%; thus 

||Q5QG(r)^(0)|U fe < \\\QSQG{t)\\\\\{Q{A + C)Q)- 1 QSVw\\h^ 

fs-f-x 

< Ke-^y \\ H 2 +i , 
provided V G H^ +2 . Concerning the second term, we obtain 



QSQ [ G(T-a)QSQG(a)ip da 
Jo 



< \\\QSQ\ 



G(t - a)QSQG(a)rp da 



H k+i 



<Po\\hI 







k+2 



< Ke-^\\v Q \\ Ht+2 , 

provided V G H^ +A . In conclusion, 

\\QSQMr)\\ Xk < Ke-^y \\ H 2 +2 , 

and analogously, 

\\VSQMr)\\ Xk < Le-^y \\ H 2 +2 , 

for some constant L > 0. Finally, it is possible to find constants v > and 
MfllVHjyfc+4) > such that 



||^ 1 + ^)(r)|| Xfc <M(||K||^ +4 )e-^|| Wo ||^ 



k + 2 



Coming back to the estimate (19.71) of r^, for any time t we have 



\n(t)\\ Xk = CMe f e- D ^\\w,\\ H 2 + ds < C|KH^ + / 



Lemma 9.2 Under the same assumptions of Proposition EH for any T, < 

T < oo, there is a constant C independent of e such that 

\\r b (t)\\ Xk < Ce 2 , (9.10) 

uniformly for < t < T . 
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Proof. The error of the bulk part of the asymptotic expansion satisfies ( 19.31) with 
/ = and, instead of g, 



9b(t) 



dt 



QSQMt) 



Since the inhomogeneous term g b (t) has a non uniform behaviour with respect 
to e for small times, we split the inhomogeneous term g(t) into the sum of two 
functions, say g b Q and g b i, as follows 

9w(t) = Ve9b(t) , = g b (t) - g b0 (t) , 

where rj e is a not increasing C°°-function such that 

f 1 for i< e/2 , 

[0 for t > 3e/2 . 

We write the error r b as the sum of two parts r b = r b o + r b i, solving the equation 

dr b0 



dt 



Sr w + -(A + C)r b0 + eg m , r b0 (0) = 
e 



and an analogous one with the inhomogeneous term g b \. Concerning the error 
r b0 , the following estimate holds by using Prop. 18.21 



r b0 (t)\\x k < Ke 



< Ke 



/3e/2 n3e/2 / 

Uo(s)\\ Xk ds < Ke f 



r 






ds {S) 



3e/2 



hollfli ds < K\\w \\ H 4 e 2 . 

fc+1 fe+1 



Finally, we consider the evolution equation for r b \. we decompose again such an 
error as 

ni = hi + h{t) , 

by introducing the auxiliary function h, which solves the problem 

— = ~Q{A + C)Qh + eg bl , h(0) = . 
Consequently, the function f fel satisfies the initial value problem 

?D± = Sr b i + -(A + C)hi+SQh, r 61 (0) = 0, 
at e 

thus it can be easily estimated in terms of the auxiliary function h as 



\W)\\x h 



< 



\SQh(s)ds\\x h ds . 
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Again by the properties of the operator (A + C), the solution reads as follows 

for t<e/2, 

e Sl/2 G ^ ~ s )9bi(s)ds for t > e/2 , 
with G e {r) bounded semigroup generated by (l/e)Q(A + C)Q. 



h(t) 



\hi(t)\\x k < e [ \\SQh(s)ds\\ Xk ds 

J e/2 

< eK f f e-»^\\g bl (s')\\ Ht ds'ds 

J e/2 J e/2 

< eK f f e ~ Uk+1 ^ (l + \\wq\\ H 4. ^ 

J e/2 J e/2 \ 

< K \\ w o\\h}, e 2 , 



s J k+1 



by applying again Prop. I8.2[ for any t G [0, T], where the constants K depend on 
T. In conclusion, 

\\r b (t)\\ Xk < K\\w \\ H 4 +i e 2 . 
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